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Entanglement is known to be an essential resource for many quantum information processes. How¬ 
ever, it is now known that some quantum features may be acheived with quantum discord, a gen¬ 
eralized measure of quantum correlation. In this paper, we study how quantum discord, or more 
specifically, the measures of entropic discord and geometric discord are affected by the influence 
of amplitude damping decoherence. We also show that a protocol deploying weak measurement 
and quantum measurement reversal can effectively protect quantum discord from amplitude damp¬ 
ing decoherence, enabling to distribute quantum correlation between two remote parties in a noisy 
environment. 


I. INTRODUCTION 

Quantum correlations are essential resources that make 
various quantum informational processes possible, and 
quantum entanglement has been in the vanguard due to 
its fundamental roles in non-locality and advantages in 
many quantum information processing mm- However, 
entanglement is not the only quantum correlation. 01- 
livier and Zurek proposed another type of quantum cor¬ 
relation, now known as quantum discord, from the per¬ 
spective of information theory [3] . Quantum discord is a 
measure of nonclassical correlations between two subsys¬ 
tems of a quantum system. The correlations arise from 
quantum physical effects. However, it does not neces¬ 
sarily require quantum entanglement. Hence, there ex¬ 
ist separable states with non-zero discord. There have 
been significant efforts made to understand the opera¬ 
tional meanings of quantum discord HIS], and find its 
applications in quantum information processing [MSI- 
Quantum correlations, both entanglement and quan¬ 
tum discord, can be degraded by decoherence which is 
often caused by unavoidable coupling with the environ¬ 
ment. There have been many studies that attempt to 
protect entanglement by tackling decoherence. For exam¬ 
ple, one can distill a highly entangled state from multiple 
copies of partially entangled states [SHU]. Decoherence- 
free subspace [131 [IH and quantum Zeno effect [T^ can 
be also used to cope with decoherence. Recently, it has 
been shown that the weak measurement and its reversal 
measurement can effectively protect entanglement from 
the amplitude damping decoherence [T5]. Many of these 
protocols might be suitable for protecting quantum dis¬ 
cord, however, no quantitative research has been done to 
show the feasibility of the protection of quantum discord. 
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Since quantum discord can exist without entanglement 
and it provides quantum advantages, protecting quan¬ 
tum discord can be useful for some quantum information 
tasks. 

In this paper we revisit the original protocol that uti¬ 
lizes the weak measurement and quantum measurement 
reversal in order to supress the effect of decoherence [T6l - 
m and investigate the protocol in terms of quantum dis¬ 
cord. We theoretically and experimentally evaluate the 
effectiveness of quantum measurement reversal in pro¬ 
tecting the amount of quantum discord. Our results ulti¬ 
mately verifies that general quantum correlations can be 
protected by the protocol. 

The remainder of the paper is organized as follows: Af¬ 
ter a brief review of quantum discord in Section H, we 
provide a numerical method to estimate quantum discord 
from a given density matrix in detail in Section HI. Then, 
we introduce the weak measurement and quantum mea¬ 
surement reversal protocol as well as the simulation re¬ 
sult on quantum discord in Section IV. The experimental 
setup and discussion is provided in Section V, and finally, 
in Section VI, we summarize our research and conclude. 


II. QUANTUM DISCORD: THE DEFINITION 

There exist variant versions of quantum discord, which 
will be introduced and discussed in the following subsec¬ 
tions. 


A. Entropic discord 

For a classical system, information entropy or the Shan¬ 
non entropy measures the ignorance about a discrete ran¬ 
dom variable X with possible values {xi,X 2 , ■■■,Xn}- If 
the probability mass function is defined as P{xi), then 
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the Shannon entropy is defined as follows [31 HD] : 

H{X) =Y^P{xi)I{x,) = -Y^P{x,)\ogf,P{x,), (1) 

i i 

where I is the information content of X, and b — 2 for 
bit. Using the definition of the Shannon entropy, we can 
find the mutual information of two random variables A 
and B, 

I{A:B)^H{A) + H{B)-H{A,B), (2) 

where H{A,B) denotes the joint entropy of two random 
variables A and B. 

The quantum equivalence of information entropy and 
mutual information are similar to their classical counter¬ 
parts. In quantum information theory, the entropy of a 
density matrix p is given by the von Neumann entropy, 

'S'(P) =-Tr(plogbp). (3) 

Note that, for a qubit, b = 2 since this normalizes the 
maximum entropic information of a qubit to 1. For a 
joint density matrix the mutual information I{pab) 
shared by quantum systems A and B is given by the 
following equation; 

KPab) = ^{Pa) + S{pb) — (4) 

where Pa{Pb) can be deduced by the partial trace 
TrB(A) Pab- In. order to get the amount of quantum 
discord D{pab), one needs to deduct the measure of cor¬ 
relation in the classical limit J{Pab) from the mutual 
quantum information I{Pab)'- 

D{Pab) = KPab) ~ J{Pab)^ (5) 

J{Pab) = sup I{PAB\{Bk}), (6) 

{Bk} 

where {Bk} is a measurement performed locally on the 
system B. 

It is noteworthy that quantum discord is not gener¬ 
ally symmetric under the exchange of the local system 
measurements. For instance, if we can perform a set of 
measurements {^fc}, instead of {Sfe}) then we may get a 
different amount of quantum discord. Note that a sym¬ 
metric discord has been proposed in order to ensure the 
symmetry [21]. Nonetheless, this paper follows the tradi¬ 
tional definition of quantum discord, because the system 
of interest generally considers the environment that has 
symmetric effects on the systems A and B. 

B. Geometric discord 

Because we need to find the supremum of I{PAB\{Bk\), 
the quantum discord between systems A and B is not 
trivial to calculate. In fact, except for special classes of 
states such as two-qubit X density matrices, there does 


not exist a closed form solution for quantum discord |221 
I23j . As a consequence, one needs to implement complex 
numerical methods in order to calculate the amount of 
quantum discord which is presented in Sec. III. 

In order to overcome this problem, Dakic et al. in¬ 
troduced geometric quantum discord that is based on 
the Hilbert-Schmidt distance between the density matrix 
Pab and its closest classical state Pab^ I-6-i B){p\b) = 
0 [Ml im • Its definition is as follows: 

Dg{Pab) = inf WPab “ PabWi^ (V 

{Bk} 

where ||X||i is the Hilbert-Schmidt 1-norm, defined as 
||A'||i = tr(VX^X). This definition of quantum discord 
also requires numerical methods. However, the calcula¬ 
tion process is much simpler and faster compared to the 
entropic definition of quantum discord since there is no 
need to perform logarithms of matrices. 

There is another definition of geometric discord, based 
on the Hilbert-Schmidt 2-norm, 

BqHpab) = inf WPab ~ PabWIj (8) 

{Bk} 

where ||X ||2 = y^tr{X^X). However, recently it has been 
pointed out that this definition is not a good measure of 
quantum correlation, because it may increase under local 
reversible operations on the unmeasured subsystem [26j . 
Hence, the discussion about the 2-norm definition will be 
omitted in this paper. 

III. NUMERICAL METHODS EOR QUANTUM 
DISCORD ESTIMATION 

In this section, we provide numerical methods to cal¬ 
culate two different definitions of quantum discord, en¬ 
tropic discord and geometric discord. For simplicity, the 
discussion starts with two-qubit density matrix (2 0 2 
systems), and we extend the discussion further for any 
multi-qudit systems {d®d’ systems). Note that the com¬ 
putational complexity of quantum discord is classified as 
NP-complete |23j . Hence, resources required for comput¬ 
ing quantum discord grow exponentially with the dimen¬ 
sion of the Hilbert space. For any d ® d' systems, we 
layout a numerical recipe for computing quantum dis¬ 
cord based on the Monte Carlo sampling of the d- and 
d'-dimensional spaces. Our method does not search over 
the entire Hilber space, but it does give us reasonably 
close results, as we tested the integrity of the algorithms 
with repetitive trials of randomly generated density ma¬ 
trices with known analytical solutions. 

A. Entropic discord 

The estimation of entropic discord consists of two parts. 
One part is to calculate I{Pab) (Ed- and it is fairly 
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trivial. Note that is in the basis of \i) ® \j) or \ij), 
where i,j S {0,1}. The other part is to find the supre- 
mum of the functional J{PAB\{^k}) (Eq. (©)■ Equation 
can be expanded to a more explicit form [3], 

JiPAs) = sup {S{pa) - S{pAB\{Bk}))- (9) 

{Bk} 


The second term in the equation is what requires a nu¬ 
merical approach. Let us define the second term of Eq. 
([^ as a function, 


HPab) 


ini S{p^B\{Bk})- 

x^ki 


( 10 ) 


We are looking for the infimum of the functional because 
Eq. (|^ has to be maximized. 

A qubit can have outcomes of either |0) or |1). How¬ 
ever, any rotational transformation of |0) or |1) is a valid 
outcome of the measurement as well. Hence, for this 
calculation, we need to consider all the possible measure¬ 
ment basis. 

We start with two orthogonal measurement bases IIo 
and Hi, 


no = 


1 0 
0 0 



( 11 ) 


By a simple rotational transformation V, we can gener¬ 
alize the measurement outcome {Bk}. 

V(6»,0) = ^(I-wl(6<,(/))cr), (12) 


Bl = vin,v, * e {0,1}. 


(13) 


Note that d is a unit vector in the Bloch sphere repre¬ 
sentation, 

( sin 6 cos (f) \ 

sin0 sinc^ j , (14) 

cos 9 j 

where 0 < 6* < tt and 0 < ^ < 27r. cr is a tensor of the 
Pauli matrices 


(15) 





j 

( CT, 




<T = 

cr2 





(^3 

where 

fo 

n 

fo 

—i 

CTi = 1 

1 

0 J 


0 


era = 


1 0 

0 -1 


.(16) 


Using the above relations, we can deduce Pab for a 
given set of measurement {5^}, 


PAB\{Bk} 


= E ( 17 ) 

• ^ rr. 11 P'i' 


*e{o,i} 


where pi is given by pi = Tr{(I (g) Bl)pj^s{I ^ Bl)}. 
It is now obvious that the functional S{pAB\{Bk}) of 
Eq. (101 is a function of 9 and (j), and we can numer¬ 
ically estimate the extremum by simply searching over 
the spherical space, 0 < 9 < tt and 0 < (j> < 2 tt. 


B. Geometric discord 


Geometric quantum discord can be calculated in a simi¬ 
lar manner. First, one needs to define an abtitrary zero 
quantum discord state for a given joint density matrix 
Pab- Eor this, let us define the reduced density matrix 
Pg given the measurement |i) oi A, i G {0,1}, i.e. Pb\\o)a 
and Pb\\i)a> 


Pb\\o)a - 
Pb||i)^ = 


1 0 
0 0 

0 0 
0 1 


(g) I 
(g) I 


Pab}) 

PabI- 


(18) 

(19) 


Then, the zero quantum discord state p^^B oan be found 
by using the following equation: 


Pab= E (vln,v)®pB|p>.- (20) 

* 6 { 0 . 1 } 


Using the relations described above, Eq. Q can also be 
calculated by searching over the same spherical space, 
0 < 9 < TT and 0 < </) < 27r. 


C. Discord estimation for arbitrary d (g) d' quantum 
systems 

Although this approach works fine, the optimization may 
be necessary for special cases. It is relatively easy to cal¬ 
culate the quantum discord of a two-qubit state, but for 
high dimensional qudits of d > 2 where d stands for a 
dimension of the quantum state, the searching process 
might take a very long time. One might encounter mul¬ 
tiple local minima for a given arbitrary density matrix. 
Though we cannot yet rigorously prove which method of 
numerical estimation is the best way to deduce the quan¬ 
tum discord of an arbitrary quantum system, a number 
of numerical evaluations led us to the conclusion that 
the Monte Carlo sampling is sufficient for estimating the 
measure. Because it may be useful for calculating the 
quantum discord for a multi-qudit system, the general 
method is briefly discussed in the following. 

For a qudit system, one can define the generalized 
Bloch sphere using the generalized Gell-Mann matrices, 
which are essentially the Pauli matrices equivalence of 
higher-dimensional extensions m- For a qudit system 
oi d = N, i.e., S\J{N), there are a total of — 1 Gell- 
Mann matrices. They can be classified into three groups: 

i) Ad-i) gynimetric Gell-Mann matrices 

= \ j){k\ + \k){j\,l <j<k<d, (21) 

ii) Ad-i) antisymmetric Gell-Mann matrices 

= -'i\j){k\ +i\k){j\,l <j<k<d, 


( 22 ) 
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iii) (d — 1) diagonal Gell-Mann matrices 


A' 



E 




|j)0'l + ^1^ + 1)(^ + i| 


,1 < Z < d-1. 


(23) 

Using the Gell-Mann matrices, we can define the gener¬ 
alized Bloch vector expansion of a density matrix 


the quantum measurement reversal [mu]- Our sys¬ 
tems of interest are two-level quantum systems (S) whose 
bases are |0)s and |l)s- Gonsidering an environment (E) 
is initially at |0)e, we can model the amplitude damping 
decoherence [T], 

|0)sG|0)£; ^|0)sG|0)e, (26) 

|l)s G |0 )e ^ Vd\1)s G |0)e + Vd\0)s G |1)b,(27) 


V =hl + Vdb-A), (24) 

a 

where the Bloch vector b = {{bi’^},{bi^},{b^j}) and A 
is the tensor of the generalized Gell-Mann matrices m- 
Let Orf be the set of all points b G such that V 

is positive semidefinite. By definition, 0 S is the 

state space or the generalized Bloch sphere. If one uses a 
systemmatic approach to calculate quantum discord as it 
is described for entropic and geometric discord for two- 
qubit states, we can search over all the generalized Bloch 
sphere, of which method consumes extensive computa¬ 
tional resources. However, for the Monte Garlo sampling, 
each component in b is just a random variable. Program¬ 
matically speaking, we select — 1 random variables 
{z/i, i/ 2 ,..., i/d 2 -i} and one additional random variable r, 
whose absolute values are all uniformly distributed in the 
range from 0 to 1. With these random variables, we can 
construct b by the following way: 


b= (25) 

By having r, we can cover all the possible Bloch vec¬ 
tor, |6| < 1. Note that these randomly chosen density 
matrices must have physical values, i.e., they must be 
Hermitian and positive semidefinite. For instance, there 
is a possibility that diagonal terms of V can be negative, 
if we carelessly applied the construction described above. 
One must be careful and eliminate such cases for calcu¬ 
lation. The estimation method of the Monte Garlo sam¬ 
pling is tested under various cases, including two-qubit, 
two-qutrit cases, and any combinations of arbitrary d- 
dimensional quantum systems. Using this method with 
a sufficiently large number of sampling can provide you 
a good estimation of quantum discord quickly. The plots 
and figures of the papers are generated with the meth¬ 
ods described in this section, including the Monte Carlo 
sampling. 

IV. THEORY 

A. Weak measurement and quantum measurement 
reversal protocol 

Let us introduce the amplitude damping decoherence 
suppression protocol using the weak measurement and 


where 0 < 13 < 1 is the magnitude of the environmen¬ 
tal decoherence and D = 1 — D. Note that amplitude¬ 
damping decoherence is a widely used model for various 
qubit systems [T] 

The experiment considers a quantum communication 
scenario depicted in the following. Alice prepares a two- 
qubit correlated state 1$), 

1$) =a|00)s + /3|ll)s, (28) 

where |Q:p-|-|/3p = 1. This state is then delivered to 
Bob and Charlie through the quantum channels of which 
amplitude-damping decoherences are characterized as Di 
and D 2 . The initially correlated state |$) is then altered 
by the amplitude-damping decoherence, and the conse¬ 
quent two-qubit quantum state shared by Bob and 
Charlie is now given as m 

f \a\^+DiD2\P\^ 0 0 y/DiDiap* \ 

^ 0 UiUzIdl^ _ 0 0 

0 0 UiUaldl^ 0 

V VDT_D2a*P 0 0 DiD2\P\^ / 

(29) 

where Dk = I — Dk, k G {1,2}. 

We can make it turn around by sequential operations 
of weak measurement (M^j^) and reversing measurement 
(Mrev), performed beforehand and afterward of decoher¬ 
ence, respectively. These operations are non-unitary and 
defined as follows : 

M„*(p.,P2) =(J 

MreviPr^^Pr^) = q 1 ) ^ ( Q 1 

where pi and Pn are the strengths of the weak mea¬ 
surement and the reversing measurement for Bob {i = 

1) and Charlie (i = 2), respectively. We chose the 
strength for the reversing measurement for protecting 
the amount of correlation of the joint state p^ to be 
Pr, = (1-A)Pz +A naiizi. Assuming that the ex¬ 
periment is performing the weak and reversing measure¬ 
ments, the two-qubit state p^ is now given as 

/ |ap-fpiP2UiU2|dP 0 0 oiP* \ 

0 PiAIdP 0 0 

a\ 0 0 P2U2|dP 0 

\ a*p 0 0 |/3|" / 

(32) 
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FIG. 1. Theoretical estimation of entropic discord as func¬ 
tions of decoherence and weak measurement; Plots (a) and 
(b) are for the maximally correlated state lil?) with |a| = |/3|, 
and plots (c) and (d) are for the non-maximally correlated 
state I'F) with |a| < \j5\ with a = 0.42. Entropic quantum 
discord under the influence of decoherence is shown in plots 
(a) and (c), whereas the effect of the weak and reversing mea¬ 
surements is shown in plots (b) and (d). Plots (b) and (d) 
are taken with Di = 0.6 and D 2 = 0.8. 


where ^ = 1 -|- {piDi{\ + P 2 D 2 ) + p2D2}\f3\'^ and pi = 
1 — Pi- Since we have the exact forms of the density ma¬ 
trices, we can analyze various quantum correlations un¬ 
der the amplitude damping decoherence with and with¬ 
out weak measurement and quantum measurement re¬ 
versal protocol. Note that the entanglement behaviour 
has been investigated in this scenario US]. The results 
showed that entanglement can be protected from the 
amplitude damping decoherence and even entanglement 
sudden death phenomenon can be avoided. 


Quantum discord protection 

We examine how entropic discord {D{p)) and geomet¬ 
ric discord {Dq{p)) behave under different decoherence, 
weak measurement, and the corresponding chosen revers¬ 
ing measurement. Note that since both and have 
forms of so called X-state, there exists an analytic solu¬ 
tion for quantum discord [22) . We have confirmed this 
analytic solution and our numerical methods in Sec. Ill, 
provide the same results. For checking the intergrity of 
our code, we have searched over tens of thousands ran¬ 
domly chosen density matrices, and it confirmed that our 
numerical method provides sufficiently close estimations, 
compared to the analytic results. Figure shows the 
entropic quantum discord and the geometric quantum 
discord, respectively. For both cases, two particular ini¬ 
tial states of |a| = |/I| and |a| = 0.42 < |/3| are inves¬ 
tigated. The plots clearly show that decoherence affects 
the two qubits independently, and their correlations can 


be circumvented by exploiting weak measurement and 
quantum measurement reversal. However, it is notewor¬ 
thy that, for quantum discord, the amplitude damping 
decoherence does not cause sudden death of correlation, 
unlike entanglement sudden death. 


V. EXPERIMENT 


Figure shows the experimental setup with photonic 
polarization qubit implementation. First, in order to 


generate two-qubit entangled state, Eq. (28) with |a| = 


1/31, type-I frequency-degenerate spontaneous parametric 
down-conversion has been implemented (not shown in 
Fig.[|). 405 nm diode laser beam is pumped into a 6- 
mm-thick /3-BaB204 crystal to generate 810 nm photon 
pairs. The down-converted photons are filtered with a set 
of interference filters whose FWHM bandwidth is 5 nm. 

There are three main parts to implement the proto¬ 
col: weak measurement, amplitude damping decoher¬ 
ence, and reversing measurement. The weak and revers¬ 
ing measurements are implemented with a set of Brewster 
angle glass plates (BPs) and half wave plates [3S|. Note 
that because the weak and reversing measurements can 
be mapped to the polarization dependent losses, it is nat¬ 
ural that the measurements can be implemented by BPs 
and half wave plates. 

The amplitude damping decoherence is implemented 
with the displaced Sagnac interferometer nil. The inte- 
ferometer couples the system’s polarization qubit to the 



FIG. 2. Experimental setup. The initial two-qubit state |<I>) 
is the two-photon polarization state. Amplitude-damping de¬ 
coherence (D) is implemented with the displaced Sagnac in¬ 
terferometers. Brewster-angle glass plates (BPs) and half¬ 
wave plates (Hs) are employed to perform weak mesurements 
(M^k) and the reversing measurements {Mrev)- Waveplates 
(WPs), polarizers (Pol.), single photon detectors (SPDs), and 
a coincidence counting unit (GGU) are used for quantum state 
tomography. 
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Decoherence, D 


Weak measurement strength, p 


FIG. 3. Experimental data for protecting entropic (top) and 
geometric (bottom) quantum discord from decoherence using 
quantum weak measurement and quantum measurement re¬ 
versal. The blue curves are for the state |a| = |/3|, whereas 
the red ones are for |q| = 0.42; (a, c) As D increases, the 
amounts of entropic and geometric quantum discord gradually 
decrease, (b, d) Even under the effect of strong decoherence 
[D = 0.6), we can reverse the amounts of quantum discords 
between Bob and Charlie by performing and Mrev(j3). 

The error bars represent the statistical error of ±1 standard 
deviation, and the dashed lines in (a) and (b) represent the 
corresponding concurrence plots. 


and quantum measurement reversal. Figure |^b), (d) 
show the entropic and geometric discords of the two- 
qubit state Pj., respectively. The reversing measurement 
parameter = p(l — D) + D is chosen for a given weak 
measurement strength p. As shown in Fig. the ex¬ 
perimental results show that the sequential operations of 
weak measurement and reversing measurement can in¬ 
deed protect quantum discord. 


VI. CONCLUSION 

We first provided numerical methods to find both en¬ 
tropic and geometric discords. By applying the methods 
to the quantum correlation protection protocol, we suc¬ 
cessfully show that quantum discord can be protected 
from decoherence by weak measurement and quantum 
measurement reversal. The protocol described in this 
paper can be applied to other types of quantum system 
beyond two-photon polarization qubits. We believe that 
this protocol is a compelling method that can be used for 
effectively handling decoherence and distilling quantum 
correlations from decohered quantum resources. 
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